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Though the problem of Hawking radiation in de Sitter space-time, in particular details 
of penetration of a quantum mechanical particle through the de Sitter horizon, has been 
examined intensively there is still some vagueness in this subject. The present paper aims 
to clarify the situation. A known algorithm for calculation of the reflection coefficient R e j 
on the background of the de Sitter space-time model is analyzed. It is shown that the 
determination of R e j requires an additional constrain on quantum numbers eR/hc » j, 
where R is a curvature radius. When taking into account this condition, the value of -Re- 
turns out to be precisely zero. 

It is shown that the basic instructive definition for the calculation of the reflection co- 
efficient in de Sitter model is grounded exclusively on the use of zero order approximation 
in the expansion of a particle wave function in a series on small parameter 1/i? 2 , and it 
demonstrated that this recipe cannot be extended on accounting for contributions of higher 
order terms. So the result R e j = which has been obtained from examining zero-order term 
persists and cannot be improved. 

It is claimed that the calculation of the reflection coefficient R e j is not required at all 
because there is no barrier in the effective potential curve on the background of the de 
Sitter space-time, the later correlate with the fact that the problem in de Sitter space reduces 
to a second order differential equation with only three singular points. However all known 
quantum mechanical problems with potentials containing one barrier reduce to a second 
order differential equation with four singular points, the equation of Heun class. 



I. INTRODUCTION 

The problem of Hawking radiation [1] , and in particular the radiation in de Sitter space-time 
[2] and details of penetration of quantum mechanical particles through the de Sitter horizon were 
examined in the literature [3-10]. Till now remains some vagueness in this point, and the present 
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paper aims to clarify the situation. 

In the paper, exact wave solutions for a particle with spin in the static coordinates of the de 
Sitter space-time model are examined in detail, and the procedure for calculating of the reflection 
coefficient R t j is analyzed. First, for scalar particle, two pairs of linearly independent solutions are 
specified explicitly: running and standing waves. A known algorithm for calculation of the reflection 
coefficient R e j on the background of the de Sitter space-time model is analyzed. It is shown that 
the determination of R e j requires an additional constrain on quantum numbers eR/hc^> j, where 
R is a curvature radius. When taking into account for this condition, the value of R e j turns out 
to be precisely zero. 

It is claimed that the calculation of the reflection coefficient R e j is not required at all because 
there is no barrier in the effective potential curve on the background of the de Sitter space-time. 

The same conclusion holds for arbitrary particles with higher spins, it was demonstrated ex- 
plicitly with the help of the exact solutions for electromagnetic and Dirac fields in 111 ]. 



The structure of the paper is as follows. In Section 2 we state the problem; some more details 
concerning the approach used could be found in 11]. 



In Section 3 we demonstrate that the basic instructive definition for the calculation of the reflec- 
tion coefficient in de Sitter model is grounded exclusively on the use of zero order approximation 
$(°)(r) in the expansion of a particle wave function in a series of the form 

*(r) = *<°)(r) + (±) *«(r) + (^Y^Hr) + ... (1) 



R?) w \R?. 

What is even more important and we will demonstrate it explicitly, this recipe cannot be 
extended on accounting for contributions of higher order terms. So the result R e j = which will 
be obtained below from examining zero-order term persists and cannot be improved. 

II. REFLECTION COEFFICIENT 

Wave equation for a spin particle (M is used instead of McR/fi, R is the curvature radius) 
reads 

' d a ^g a P dfs + 2 + M 2 ) V(x) = , (2) 



f-9 

and is considered in static coordinates 

dS 2 = <&dt 2 - — -r 2 (d9 2 + sin 2 6d(/) 2 ) , < r < 1 , $ = 1 - r 2 . (3) 
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For spherical solutions = e let f(r)Yj m (9,(/)), e = ER/hc, the differential equation for f(r) is 

+ (r + ¥) ^ + ($2 i ^r 2 " ) ^ = ° ' (4) 

All solutions are constructed in terms of hypergeometric functions (let r 2 = z) : 
regular at r = standing waves are given as 

f(z) = zi' 2 (1 - z)^' 2 F(a b,c;z ) , k = j/2 , a = -it/2 c = j + 3/2 , 
_ 3/2 + j + iy/M 2 - 1/4 -ie u _ 3/2 + j - i^/M' 2 - 1/4 - ie 

a - g ' b ~ 2 ; { ' 

singular at r = standing waves are 

g (z) = z-d+i)/ 2 (1 - z) -^ 2 F(a, 0, T ,z), k = -(j + l) /2 a = - ie/2 , 

1/2 -j + j^M 2 - 1/4 -ze 1/2 - j - j^M 2 - 1/4 - ze 
c = -J + 1/2 , a = , /3 = - . (6) 

With the use of the Kummer's relations, one can expand the standing waves into linear combi- 
nations of the running waves 

T(c)T(c-a-b) t T(c)T(a + b-c) 

J[ ) T(c-a)T(c-b) run[ ) + T(a)T(b) run{ ) ' l J 

UZ? n (z) = z?' 2 (1 - z)~^ 2 F(a, b, a + b - c + 1; 1 - z) , 
tC„ W = ^' /2 (1 " ^) + " /2 ^(c - a, c - b, c - a - b + 1; 1 - z) , (8) 



a* = ( c - a) , b* = (c- b) , (a + b - c)* = -(a + 6 - c) , 



[ c/; n t (*) ]* = tCnW , (9) 



„ , . „ „ T(c)r(c — a — 6) . . „ _ T(c)T(a + b — c) TT , . 

" 2 Re ric - W-t c, »"' (2) = 2 Re r(W) c, - <2) ' (10) 



Similarly for g(z) 



r(7)r(7 _ a _m r( 7 )r(a + /3- 7 ) 



= i> „)iv 1] «W + r(a)r(/3) • < u > 

£ClW = ^' /2 (1 - ^)"" /2 F(a + l- 7 ,/3 + l-7,a + /3 + l-7;l-^), 

CCn W = ^' /2 (1 " z) +le/2 F(l - a, 1 - /?, 7 + 1 - a - /?; 1 - z) , (12) 
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5 (,) = 2RePlfc^^) = 2Re 



r( 7 )r(a + /3 - 7 ) 



r( 7 -a)r( 7 -/3) > r(a)r(/3) 

Asymptotic behavior of the running waves is given by the relations 



U in (z) 



UZ'nir ~ 0) ~ -L , Eft*(r ~ 1) ~ (1 - r 2 )-/ 2 , 



CCn(r~0)~ 



1 



r j+i 



CC»(r ~ 1) ~ (1 - r 2 )+-/ 2 , 



or in new radial variable r* G [0, 00): 



R , 1 + r 

— In , 

2 1-r ' 



exp(2r*/i?) - 1 



" exp(2r*/i?) + 1 ' 
CC*(r* ~ 00) ~ (2-^/^ c ) exp(+iEr*/hc) , 

U™ n {r* ~ 00) ~ (2+^/^) exp(-iEr*/hc) , e = £i?/ftc 



For the standing waves we have 

/(r ~ 1) ~ 2 Re 
g(r ~ 1) ~ 2 Re 



f(r ~ 0) ~ r J ' , ^(r ~ 0) 



r j+i ' 



I^±^2^/-ex P (-^7M 



(13) 



(14) 



(15) 



(16) 



On can perform the transition to the limit of the flat space-time in accordance with the following 
rules: 



p + 1 - ieR + i^WM 1 - 1/4 



t _,, + i-fcit-y«w-i/4 



x „ 9 , s ab z a(a + l)b(b + 1) z 2 . . 

lim (i? 2 z) = i? 2 , F(a,6,c;z) = l + -- + 1 ^- >- - + ..., (18) 
i?^oo c 1! c(c + 1) 2! 



and further (i? — >■ 00) 



a + n 1 lp + 1 I ~ 1 \ n —ie — iM 

1 ie + i\ M 2 - ——- + — - 



i? 2\R V 4R 2 R 2 



6 + n 1/p+l /„,„ 1 \ n -ie + iM , N 

— =2(V- £ -V M, -4ffj + 5' S — 2— • (19) 

Let e 2 — M 2 = A; 2 , then we arrive at 

lim F(a, b,c;z)=T(l+p)Y / rn g ^ , , 
i?-^oo ^ n!r(l + n + p) 

- c + 1, „ - c + 1, -c + 2; „) = r(l - p) £ „; r(1 + ,;_' p) ■ (20) 
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Allowing for the know expansion for Bessel functions 



J p (x) 



(?) P E 



(ix/2) 



2/i 



n!r(l + n + p) ' 



we get (wave amplitude A will be determined below) 



1 



lim AU%*(z)=hm A -j=, 



T(-ieR + l)r(-p) T(l + p){2/k)PR- p+1 / 2 



T[U+i^/R 2 M' 2 - 1/4 -ieR + p + l)]T[U-i^ R 2 M' Z - 1/4 - ieR + p + 1)] 



J p (kr) 



+ 



T(-ieR + l)r(+p)r(l - p){2/k)- p R +p+l l 2 



T[U+iy/WW - 1/4 -ieR-p + l)]T[U-i^R 2 M 2 - 1/4 -ieR-p+ 1)] 



J- P (kr) 



Performing limiting procedure (see the detail in [11]) we derive the relation 



lim A Uout(z) 

R—too 



where 



H 



(i) 



(*) 



1 ,fl H W f kr ) 



e** J p (s) - J_ p (x) 



(21) 



(22) 



stands for Hankel spherical functions. 

In connection with the limiting procedure, let us pose a question: when the relation (|22p gives 
us with a good approximation provided the curvature radius R is finite. 

This point is important, because when calculating the reflection coefficient in the de Sitter space 
just this approximation ([22]) was used [3-6]. 

To clarify this point, let us compare the radial equation in Minkowski model 



d 2 



2 d 



dr 2 ^ r dr 



M 



and the appropriate equation in de Sitter model 



j(j + 1) 



M 2 + 2 



f 







j{j + 1) 



d 2 2(l-2r 2 /i? 2 ) d 

d^ 2 + r(l-r 2 /R 2 ) Ik + ~l-r 2 /R 2 ) 2 l-r 2 /R 2 r 2 (l-r 2 /R 2 ) 



At the region far from the horizon r R, the last equation reduces to 



d 2 



2 d 



dR 2 + RdR + 6 



i(j + 1) + 2 + 







(23) 



fej = 0- (24) 



(25) 



R 2 R 2 

So, we immediately conclude that eq. (f2"3"j) coincides with ([23]) only for solutions with quantum 
numbers obeying the following restriction 

2 



M 2 > 4> . 



(26) 
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In usual units, this inequality reads 

t? 2 x h 

h, = a mc , A = — 

mc 

Instead, in massless case we have 

R 2 u 2 



R 2 



or 



10" 



A 2 

- 1 » ^ f 



RHtT 2 2 



(27) 



(28) 



So we can state that the above relation (|22p is a good approximation at a finite R only for 
quantum numbers obeying ([2T|) - (pg]) . 

One additional point should be emphasized, the radial equation in de Sitter space can be 
transformed to the form of the Schrodinger like equation with an effective barrierless potential. 
Indeed, in the variable r* eq. (|24"|) reduces to 

~ d 2 



dr 



*2 



U(r*) 



G(r* 







U(r* 



1 



„2 r 



R 2 



4(1 -r) 



1 + r 



+ m 2 R 2 + 



(29) 



It is easily verified that this potential corresponds to attractive force in all space 



F r 



dU 

dr* 



1 — r 



+ 



2r 



^±11 + m 2 R 2 + 



+ 4(l-r)) + (1 



2^ /2j(j + 1) 



+ 4 



1 



1 + r 



1 + r 
> 



At the horizon, r* — > oo, the potential function U(r*) tends to zero, so G(r*) ~ exp(±ier*). 

The form of the effective Shcrodinger equation modeling a particle in the de Sitter space indicates 
that the problem of calculation of the reflection coefficient in the system should not be even stated. 
However, in a number of publications such a problem has been treated and solved. So we should 
reconsider these calculations and results obtained. Significant steps of our approach are given 
below: 

The existing in literature calculations of non-zero reflection coefficients R e j were based on the 
usage of the approximate formula for U°^(R) for the region far from horizon (|2ip . However, as 
noted above, the formula used is a good approximation only for solutions specified by (|26p . that is 
when j <C eR. 

It can be shown that the formula existing in the literature gives a trivial result when taking into 
account this restriction. The scheme of calculation (more detail see in [11]) is described below. 

The first step consists in the use of the asymptotic formula for the Bessel functions: when 
x S> v 2 we have 



J(x) 



r(2»/ + i) 2^- 1 / 2 j_ 

T(u + 1) r(i/ + l/2) 7^ 



7T 



1 



exp +i(x - -{v + -)) + exp -i(x -)) 



7T 



1, 
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so when j < j 2 <C eR « eR we derive 



' e +ieR 

eR 

ieR 



d exp(-t~(p+-)) + C 2 exp(-i|(-p+-)) 



+ 



ri? Ci exp(+*|(p+^)) + C 2 exp(+t|(-p+i)) 



where Ci and C2 are given by 



Ci 



Ci 



_ r(a + ft + i-c)r(i-c) 2-J- 1 r(2p + 1) 

~ r(6-c+l) r(o-c+l) (eR)i r(p + l/2) ' 
r(a + 6 + 1 - c) r(c - 1) (eR) j+1 r(-2p + 1) 



(30) 



(31) 



r(a) r(6) 2-j r(p + 1/2) ' 

The reflection coefficient R e j is determined by the coefficients at e~ leR /eR and e +l<E ^/ei?. It is the 
matter of simple calculation to verify that when eR 3> j, the coefficient R e j is precisely zero 



eR » i , 



i? ej = . 



(32) 



This conclusion is consistent with the analysis performed above. 



III. SERIES EXPANSION ON A PARAMETER R~ 2 OF THE EXACT SOLUTIONS 
AND CALCULATION OF THE REFLECTION COEFFICIENT 

In Section 3 we demonstrate that the basic instructive definition for the calculation of the 
reflection coefficient in de Sitter model, being grounded exclusively on the use of zero order ap- 
proximation <&(°)(r) in the expansion of a particle wave function in a series of the form (pQ), cannot 
be extended on accounting for contributions of higher order terms. 

Let us start with the solution, the wave running to the horizon 

U out (z) = T(a + b - c + 1) [ a F(z) + p G(z) } , 
a = Eil^ 3 r(C " 1} (33) 

r(6- c + i)r(o-c + i) ' p r(o)r(6) ' 1 ; 



where 



and 



F(z) = «(p-V2)/2( 1 _ z )-^/ 2 F(a, b,c;z) , c = j + 3/2 = 1 + p , 

1 + p — ie + iJm 1 — 1/4 1 + p — ie — iJm 1 — 1/4 
a = o ' 5 = o ' ( 34 ) 



G(z) = z (-P-l/2)/2 (1 _ z )+«/2 F ( a _ c + 1, & _ c + 1, 2 - c; z) , 2 - C = 1 - p , 

1 -p - ie + iy/m 2 - 1/4 1 - p - ie - iy/m 2 - 1/4 
a — c+l = , o — c+l = . (35) 
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For the following we need the expressions for all quantities in usual units. It is convenient to 
change slightly the designation: now R stands for the curvature radius 

- 2 ER R „ McR R . h 



z = w , 6= — =,,-, E = »Mc, 



m 



n 



x 



Now, the relations (|33j) - (135j) read 



-i[iR/2\ 



II R R 2 1 

a= 2 \ l+p - i n + iv-A 



F(a,b,c; 



R? 



c = 1 +p , 



R 



lR 2 1 



(37) 



and 



v I r 2 



-ifiR/2\ 



F(a-c+l,6-c+l,2-c; -=■) 



2-c=l-p, 



'i? 2 1 



o-c+i = i[i- P -i^+iM-i), b-c+i= i - [i- P -i^-i^-z 



(38) 



The task consists in obtaining the approximate expressions for F(r) and G(r) in the region far 
from horizon, r R; first let us preserve the leading and next to leading terms (we have a natural 
small parameter X/R). 

First, let us consider the exponential factors 



2 \ ±i)iR/2\ 

=exp 



±« — r hi 1 — — „ 
2X K R 2> 



cos 



^A ln(1 "^ 



± i sin 



^A ln(1 -^ 



(39) 



Using the expansion for the logarithmic function 

x 2 x 3 . , . r 2 , / r 2 1 r 4 1 r 6 

ln(l-x) = -(, + T + y + ...), H l--) = -l— + -^ + -- + 



we get (assuming that r 2 <C Ai?) 





lr 4 1 r 6 \ 
^R 2 + 2"R4 + 3fl6 + 

r 2 1 r 4 1 r 6 \ / /u 2 r 4 \ . pR f r 2 1 r 4 \ 

fi 2 + 2"^ + 3^6 +- ) w I 1_ SA 2 ^ 2 J T% 2A I ~R2 + 2^R4 ) • 1 lm 
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Terms in eq. (|40p can be written in descending order 



±i[iR/2\ 



1 2 4 

r u r r 
1 + l i l — n -F -A 



2Ai? 2 4A 2 i? 2 



4A 2 i? 2 



A 



2Ai? 



< 1 



(41) 



Now we turn to the hypergeometric function F(a,b,c; z) from (|37p . Because i? ~ 10 30 ,A ~ 
10~ 12 , one can use the approximation of a leading and two next order terms in the expressions for 
the following parameters 



1/ R R X 2 \ 1 + p n-\R A 

a = - 1 + V — ILL— +1-1 1 —77 = I I — , 

2 I A A V 4i? 2 2 2 A 16R 



6- 2 1+P *M A 4i? 2 " 2 * 2 A+'iei? 



(42) 



Then, the hypergeometric function is given as 



„, , r 2 . 1 06 o 1 1 0(0 + 1)6(6 + 1). 2n2 

^' ' C ' F> = 1 + ^ + 2itf c(c + l) (r ) 

, 1 1 a(a + l)(a + 2)6(6 + 1)(6 + 2) 2 3 

— (r r + ••• 



SIR 6 



c(c+ l)(c + 2) 



| 1 1 a(a + l)(a + 2)...(a + ra- 1)6(6 + 1)(6 + 2)..(6 + n - 1) 2 



(43) 



n!.R 2 ™ c(c + l)(c + 2)...(c + n- 1) 

It is convenient to introduce a shortening notation for small quantity X = X/R), then a typical 



term is represented as 



1 (a + n)(6 + ra) 



1 



2A 



8(0-1) 



R 2 (c + n) p + 1 + n 

-*(A* + 1) / | , l+p + 2n ^ 



X 2 



2A 



+ 1 



8(0 + 1) 



1 



, t pdzi) 

p + l + n \ 4A 2 y 
Below we will use the notation 



2m , . (1 +p + 2n) 2 - 1/4 2 
1 + -^-^(l+p + 2n) X- V ^ „ 4 — X 2 



2 - 1 



(44) 



■2_M 2 -1 
" " A 2 ' 



then 

1 (a + ra)(6 + n) 
i? 2 (c + n) 



p + 1 + n \ 4 



2iu , , (1 + p + 2n) 2 - 1/4 9 
l + ^r-^(l+p + 2n) X- V ^ „ 4 — X 2 



2 -l 



2 -l 
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Thus, we have the following approximate expressions for the first few terms in the series 
r 2 ab if k 2 r 2 



R 2 



p + 1 



r 2 (o + 1)(6+ 1) ^ 1 

r 2 (F+T) 

r 2 (a + 2)(b + 2) 
R 2 (c + 2) 

r 2 (q + 3)(6 + 3) 
i? 2 (c + 3) 

r 2 (a + 4)(6 + 4) 
^2 ( c + 4 ) 



P + 2 
1 

P + 3 
1 

p + 4 
1 

p + 5 



4 

/c 2 r 2 

fc 2 r 2 ' 

k 2 r 2 

/c 2 r 2 



1 + ^ (1+rtA -_(l±ifcIZl x2 



1 + 



/i 2 - 1 
2ifi 



// 2 -l 



/7/ 



1 



1 + 



1 + 



1 + 



2i/j, 

pT^i 

2i[i 



(l+p + 2 x 1)X- 
(l + p + 2x2)Jf- 
(l+j) + 2x3)I- 
(l+p + 2x4)X- 



(1+P + 2X l) 2 -l/4 x2 



^ 2 - 


1 




(1 + p + 2 x 


2) 2 


-1/4 




1 




(1 +p+2x 


3) 2 


-1/4 


p 2 - 


1 




(1 + p + 2 x 


4)2 


-1/4 




1 





X 2 



X 1 



x z 



r 2 (a + n)(b + n) 
~R? (c + n) 



1 



k 2 r 2 



p + 1 + n 



1 + J^ (1 + p + 2n) x-" + '' + 2 "> 2 - 1 / 4 ^ 



//- 1 ' /i 2 - 1 

Now we are ready to write down the expressions for the coefficients of the hyper geometric series 
preserving only leading and two next order terms 



1 r 2 ab ^ 1 
V.R 2 ^ ~ p+l 



k 2 r 2 



1 + 



2ip 



{l+p)X 



(l+p) 2 -l/4 
li 2 - 1 



1 (r 2 ) 2 a6(a + 1)(6 + 1) 
2! (W c(c+l) 



/c 2 r 2 \" 1 f 2i[i 



2 



2!(p + l)(p + 2) 



I 1 + ^3y[( 1+ ^) + (1 +P + 2 x 1)] X- 



-X 2 



4 ^ r (l + p)(l +P + 2xl) + ( 1+ f -V4 + (l + P + 2xl) 2 -l/4 



( /U 2_ 1 )2V /VV ^ / /i2 _ 1 jU 2_ 1 



1 (r 2 ) 3 a6(a + l)(6+l)(a + 2)(6 + 2) / fc 2 r 2 \ 3 



3!(i? 2 ) 3 c(c+l)(c + 2) V 4 / 3!(p+l)(p + 2)(p + 3) 

x (l + I^-[(l+p) + (l+p + 2x 1) + (l+p + 2 x 2)]- 



-X 



2 V 



(m 2 - 1) 2 



[(l+p) + (l+p + 2x l)](l+p + 2x2)+ 



(1 + p) 2 - 1/4 (l+p + 2 x l) 2 - 1/4 (l+p + 2 x 2) 2 - 1/4 
fj, 2 — 1 fj, 2 — 1 /i 2 — 1 
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1 (r 2 ) 4 ab(a + !)(& + l)(a + 2)(b + 2)(a + 3)(fe + 3) 
4!(i? 2 ) 4 c(c + l)(c + 2)(c + 3) 

/ fcvV i 



^ 4 y 4!(p + l)(p + 2)(p + 3)(p + 4) 
+ X -^t[(1+p) + (1+P + 2 x l) + (l+p + 2x2) + (l+p + 2x 3)]- 



2 V 



- 1)2 



[(1 +p) + (1 + p + 2 x 1) + (1 + p + 2 x 2)](1 +p + 2x3)+ 



(1 +p) 2 - 1/4 (1 + p + 2 x l) 2 - 1/4 (1 + p + 2 x 2) 2 - 1/4 (1 + p + 2 x 3) 2 - 1/4 
H 2 — 1 fi 2 — 1 fi 2 — 1 fi 2 — 1 



1 (r 2 ) n o6(a + l)(6 + l)....(a + n-l)(6 + n-l) / fcV 



n 



! (i? 2 )™ c(c+ l)...(c + n) 14/ n!(p + l)(p + 2)...(p + n) 



-X 



2i[i 

4^ 2 



x |l + I T ^[(l+p) + (l+p + 2xl) + (l+p + 2x2) + ... + (l+p + 2x(n-l))]- 



(/x 2 - l) 2 



[(l+p) + (l+P + 2xl) + ... + (l+p + 2x (n-2))] (l+p + 2x (n-l))+ 



(l+p) 2 -l/4 {l+p + 2 x l) 2 - 1/4 (1 +J9 + 2 x (n- l)) 2 - 1/4 

fi 2 — 1 fi 2 — 1 ^i 2 — 1 

Thus, initial exact hyper geometric function can be approximated by the sum of three series 

2 

F = F(a,b,c; ^) = F (r) + X F 1 (r) + X 2 F 2 (r) . (45) 

The leading series Fo(r), in fact, reduces to the Bessel function 

F(r) = l+ {lkr/2)2 + {ikr,2)A + + (*W2) 2 " 
01 ' n!(p+l) 2!(p+l)(p + 2) "' n!(p+l)(p + 2)...(p + n) 

= r(P + i)| J' ( ^f;„ r r(i +p ) (£)-\(*0 . (40 



where 



JJx) = (-YY ,l/ 2/4) " : , x = A:r. 
pv 7 v 2 7 Y n! r (^ + 1 + n ) 

The second series is given by 



that is 



I , ,2 2 , ,x 3 [(1 + P) + (1 + P + 2 x 1) + (1 + p + 2 X 2)] 
1 7 J 3!(p+l)(p + 2)(p + 3) 

, r _, 2 2,m [(1 + p) + (1+P + 2x l) + (l + p + 2x2) + (l+p + 2x3)] 
1 /j 4!(p + l)(p + 2)(p + 3)(p + 4) 

/ /;2r2/ir [(l+p) + (l+P + 2x !) + ... + + 2 x (n -!))] ) 
1 7 J n!(p + l)(p + 2)...(p + n) J ' 



l) + (l + p + 2x2)] | 



2i(j, 


/-fcV 


/i 2 -l 


V 4 




-k 2 r 2 V 









3\(p + l)(p + 2){p + 3) 



f -k 2 r 2 \ 3 [(l+p) + (l+p + 2x l) + (l+p + 2x 2) + (l+p + 2x3)] 
+ l 4 j 4!(p+l)(p + 2)(p + 3)(p + 4) + "' 



+ 



-k 2 r 2 \ n [(1 +p) + (1 +p+ 2 x 1) + ... + (1 +p + 2 x n)] 



(n + l)!(p + l)(p + 2)...(p + n + 1) 



or shorter 



yin „ Mji „, ^,J-k 2 <- 2 \ ^. f-kVy [(l+p) + ... + (l+p + 2xn)] 
XFl{r)=X — r(p+1) ^—j E(—J („+!)! r (p + 2 + n) 



Using known sums 

[(1 +p) + ... + (1 +p + 2 x n) = (1 +p)(n + 1) + 2(1 + 2 + 3 + ...n) 

= (p + l)(l + n ) + 2^ — = (n + l)(n + l+p) 



we derive 



/i 2 -l V 4 / (ra + 1)! T(p + 2 + n) 



/i 2 - 1 v 4 y ^ n!r(p+l+n) 

/_ k 2 2\ 2 - /kr\- p , . „{-k 2 r 2 \ 2iu -,_ 



Thus, the approximation (J45J) can be presented as follows 

F = F (r) + XF x (r) + X 2 F 2 (r) = F (r) + X (=^) F (r) + X 2 F 2 (r) 

[i z — 1 4 

Similar relations can be derived for the hypergeometric series G(r): 



2 

G = F(a - c + 1, b - c + 1, 2 - c; ^ ) = G (r) + X Gi(r) + X 2 G 2 (r) . 
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The leading term again reduces to the Bessel function 



' I~ r \+P 

Go(r) = T(l-p) ( — \ J- P (kr). 



(53) 



Next order term is given by 



k 2 r 2 \ 2i(j, 



XG x (r) =X 



4 / ix l - 1 



Go(r) 



(54) 



So, the approximation (|52|) is presented as 



G = G (r) + XGi(r) + X 2 G 2 (r) = G (r) + X (-^) G (r) + X 2 G 2 (r) , 

// z — 1 4 



(55) 



Now let us consider the whole function 



-ifiR/2\ 



F(a,b,c; —5) 



R - P+ i/2 r 
\/r 



2 2 4 A 



2Ai? 2 4A 2 i? 2 ^ 4A 2 i? 2 



F (r) + 



2i/i —A; 2 ? 



/i 2 - 1 v 4 



XF (r)+X 2 F 2 (r) 



*b(r) + (-^) XFo(r) + X 2 F 2 (r)+ 

ix^ — l 4 



r 2 - r 2 2i^ —k 2 r 2 



)XFo(r)+^ — X 2 F 2 (r)- 



+^ * 4A2fi2 F o(0 + »> X 



r 2i^ —kr 



AX 2 R 2 li 2 - 1 v 4 



(— — ) **b(r) + V * 2 F 2 (r) 



r -p+i/2 J_ 
\/r 
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X 



\3Su i ,■„/■ ' \2 V 3 



2\R' uv ' ^ y 2XR 

Preserving only first two terms we have 

, r +p 

F(z) = R-P+ 1 ' 2 — 
y/r 

Similarly, for G(r) we obtain 

G{z) = R+P+ 1 / 2 T-L G (r) + V 



j F (r) + ifi 



2XR' 



X a F 2 (r) 



2„2 



A 



\ 2 R? G ° (r) + i^ G2(r) 



(56) 



(57) 



One should emphasize one feature of the expansions (|56j) and (|57|) : these approximations are real 
valued as we must expect remembering on relations (|10p and Q13p. 

In the known method of determining and calculating the reflection coefficients in de Sitter 
model [3-6], authors used the only leading terms in approximations (|56p and (I57j) . because only 
these terms allows to separate elementary solutions of the form e ±ikr - see (|30p . 

Let us perform some additional calculations to clarify the problem. The functions F(z) and 
F{z) enter the expression for (to horizon) running wave 



U out (z) = r(o + b - c + 1)[ a F{z) + (3 G(z) } , 



a 



r(i-c) 



r(6-c+l)r(a-c+l) ' 
where (introducing a very large parameter Y = R/2X) 

r(-p) 



r(c-i) 
r(a)r(6) ' 



(58) 



a 



/3 



r[-t(Ai - 1)^ + (1 - P)in V[-i{pL + 1)Y + (1 - p)/2)] 



r(+p) 



(59) 



r[-i(M - i)y + (l + P )/2)] rH(/x + i)y + (i + P )/2)\ 

For all physically reasonable values of quantum numbers j (not very high ones) and values of \i 
(different from the critical value \x = 1 and not too high ones - they correspond in fact to energies 
of a particle in units of the rest energy) the argument of T-functions in (|59p are complex-valued 
with very large imaginary parts. 

Let us multiply the given solution ([58]) by a special factor A which permits us to distinguish 
small and large parts in this expansion: 

T(a - p/2 + l/4)r(6 - p/2 + 1/4) 



.4 



T(a + b-c+ 1) 



(60) 
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So, instead of flooj) . we get 



a = T(l - c 



/3' = r(c-l) 
Instead of ([59]) . the expressions for a' and /3' are 

r[-*( M -i)r + i/4] 



AU out (z) = a' F(z) + G(z) , 
T(q - p/2 + l/4)r(b - p/2 + 1/4) 

r(6-c+i)r(o-c+i) 

r(a-p/2 + l/4)r(6-p/ + l/4) 



r(a)r(6) 



a « r(-p) 
/?' ~ r(+ P ) 

Allowing for identities 

r(-p) 



r[-i(M+ l)F + l/4] 



r[-i(// - i)y + (i - p)/2)] r[-i(ji + i)Y + (i - p)/2)] 
r[-i(/j-i)y + i/4] r[-i(/i + i)y + i/4] 



rH( M - \)y + (i + p)/2)] r[-iOx + i)Y + (i + p)/2)] 



7T 



1 



sinp7r T(l + p) ' 
a' and /3' (|62p are transformed into 

i r[-*(M-i)r + i/4] 



r(p) = + 



7T 



1 



sinp7r r(l — p) 



7T 



T[-i(n + i)y + 1/4] 



sinpTr T(l + p) T[-i(n - 1)Y + (1 - p)/2)] r[-i(M + 1)Y + (1 - p)/2)] 

7T i r[-i(p - i)y + 1/4] r[-i( M + i)r + i/4] 



sinpvr r(i - p) r[-t(// - i)y + (1 + p)/2)] r[-i(jj, + i)y + (1 + p)/2)] 

Now we have to take into account the asymptotic formula for T-function 



r(z + A) = z a~b + l (A-B)(A + B + 1) 
T(z + B) \ z 2 



then (remembering that Y = R/2X) 

rH(/i-i)y + i/4] 
r[-i( M -i)y + (i- P )/2)] 

rH(/x + i)y + 1/4] 



— i 




P/2-1/4 r 




2A ^ 





r[-i(M + i)y + (i- P )/2)] 



2A 



p/2-1/4 



and 



r[-4-i)y + i/4] 
r[-f(/i-i)y + (i+p)/2)] 

rH(/i + i)y + 1/4] 
rH(/. + i)y + (i + p)/2)] 





-p/2- 


-1/4 r 


[ 2X . 






1. 2 2A * 


-p/2- 


-1/4 , 









1 + 



1 + 



arg z \ < 7T , 


| z — )■ oo 


2A (2p - 


l)(7-2p)l 


-i{n-l)R 


32 


2A (2p - 


l)(7-2p)l 




32 


2A (-2p 


-l)(7 + 2p) 


-i{n-l)R 


32 


2A (-2p 


-l)(7 + 2p) 



-i(ji + l)R 



32 



(61) 



(62) 



(63) 



(64) 



oo ; (65) 



(66) 



.(67) 



Substituting (|66|) and ([67j) into (|Mj) . we obtain 



a: 



M 3 -l 



p/2-1/4 



i? 2 



1 + 



V 4A2 

2A (2p-l)(7-2p) 



(_l)P/2-l/4 



7T 



1 



-i(fi-l)R 



32 



1 + 



sinp7r T(l + p) 
2A (2p-l)(7-2p) 



/ 2 r \ -p/2-1/4 

^ + (^1^ ,_-n-p/2-i/4 - 



X 



1 + 



2A 



V 4A2 

-2p-l)(7 + 2p) 



-i(/U+l)i? 
(-1)" 



32 



1 



1 + 



2A 



sinp7r r(l — p) 
-2p-l)(7 + 2p) 



-i(n-l)R 32 
Preserving the only terms of first two orders we have 

7T i m+p- 1 / 2 



sinp7r T(l + p) \2 



i? +p-l/2 ( _ 1)P /2-l/4 j + 



+ 32 

(2p-l)(7-2p) /i A 



7T 



sinp7r r(l — p) \2 



8 /i 2 - 1 i? J ' 

i? -p-i/2(_ 1 ) -p/2-1/4 A _ , (~2p-l)(7 + 2p) /x A 
\ 8 [I 2 — 1 R 



Substituting expressions ([69]) into the following expansion 

AU ou \z) = a' F{z) + p' G{z) , 
F{r) = R-P^/ 2 T{l+p){^j P ^J P (kr) 

G{r)^W +l l 2 r(1 _ p) Qy_Lj_ p(A;r) 
we arrive at the following zero-order approximation 

AU ou \z) = $J"*(r) = 

1 /£A +p - 1/2 



7T 



+ 



sinp-7r r(l + p) \2 

i m-"- 1 / 2 



7T 



sinp7r r(l — p) \2 



i?-'- 1 2 ( - J )''■'' 2 ~ 1 B~ l ' +l/2 \'i\ ! (|) P -Lj p ( k r)+ 



k\ p 1 



that is 



out / 



7T /2fc 



sinp7r V r 



,(_l)P/2-l/4 Jp(fcr) + (.^-P/2-1/4 J_ p ( fcr ) 



■(-^f7 W-J_ p (At)] ; 
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which coincides with a spherical wave propagating in Minkowski space from the origin, expressed 
through the Hankel functions of the first kind 

= [{-IfUkr) - J- P (kr)} . (72) 

Sill JJ7V 

IV. CONCLUSION 

The last but not the least mathematical remark should be given. All known quantum mechanical 
problems with potentials containing one barrier reduce to a second order differential equation with 
four singular points, the equation of Heun class. In particular, the most popular cosmological 
problem of that type is a particle in the Schwarzschild space-time background and it reduces to the 
Heun differential equation. Quantum mechanical problems of tunneling type are never linked to 
differential equation of hyper geometric type, equation with three singular points; but in the case 
of de Sitter model the wave equations for different fields, of spin 0, 1/2, and 1, after separation 
of variables are reduced to the second order differential equation with three singular points, and 
there exists no ground to search in these systems problems of tunneling class. 
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